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Abstract 

Optical spectropolarimeters can be used to produce maps of the surface magnetic fields of stars 
and hence to determine how stellar magnetic fields vary with stellar mass, rotation rate, and 
evolutionary stage. In particular, we now can map the surface magnetic fields of forming solar-like 
stars, which are still contracting under gravity and are surrounded by a disk of gas and dust. Their 
large scale magnetic fields are almost dipolar on some stars, and there is evidence for many higher 
order multipole field components on other stars. The availability of new data has renewed interest 
in incorporating multipolar magnetic fields into models of stellar magnetospheres. I describe the 
basic properties of axial multipoles of arbitrary degree ^ and derive the equation of the field lines in 
spherical coordinates. The spherical magnetic field components that describe the global stellar field 
topology are obtained analytically assuming that currents can be neglected in the region exterior 
to the star, and interior to some fixed spherical equipotential surface. The field components follow 
from the solution of Laplace's equation for the magnetostatic potential. 
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I. INTRODUCTION 



The solution of Laplace's equation by separation of variables for the electrostatic po- 
tential in a region external to a charge distribution is a standard topic in graduate and 
undergraduate courses in electromagnetism.'^ As an example application, the long-range 
interaction between molecular charge clouds can be determined via a multipole expansion 
of the electrostatic potential obtained by solving Laplace's equation.^ The solution for the 
magnetostatic potential in a region devoid of current sources is equivalent to the electrostatic 
caseP 

In this paper I show how this approach can be adapted to construct models of the large 
scale magnetic fields of stars and planets. I will derive the equation for the field lines of an 
axial magnetic multipole of arbitrary degree I with a source surface. The source surface is 
a spherical surface of radius Rs used by the solar physics community to mimic the effects of 
the solar wind. The solar wind is continuous streams of outflowing charged particles that 
open the large scale magnetosphere (the region external to a star, or equivalently a planet, 
consisting of closed magnetic field lines) of the Sun.^- At Rs the field is assumed to be 
purely radial. This model has been successfully adapted to produce models of stellar mag- 
netospheres via field extrapolation from observationally derived magnetic surface maps.'^'^'^ 
As shown in Fig. [T] multipole magnetic fields with a source surface boundary condition in- 
corporate regions of closed field line loops, as well as regions of open field lines along which 
outflows (stellar winds) are launched. A basic assumption of the model is that currents can 
be neglected in the region external to the star, and interior to Rs. With this assumption and 
with the source surface boundary condition, the complicated problem of solving Poisson's 
equation to determine the magnetic field components is bypassed. Instead, the B-field is 
derived from the general solution of the familiar Laplace's equation. 

The layout of the paper is as follows. In Sec. [IT] some recent results from the study of 
stellar magnetic fields are presented. I define the field components of a multipolar large scale 



stellar (or equivalently a planetary) magnetosphere in Sec. Ill In Sec. Ill A I discuss how 
models of stellar magnetospheres can account for the distortion of the large scale magnetic 
field caused by stellar outfiows (winds) by the use of the source surface boundary condition. 



and demonstrate in Sec. Ill B how the field components are modified. The modified field 
components are derived from the solution of Laplace's equation in spherical coordinates. 
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subject to the assumptions discussed at the beginning of Sec. |IIIB[ In Sec. [IV] I solve the 
differential equation for the path of the field lines for an arbitrary multipole I with a source 
surface. I summarize the main results in Sec. IVl 

Throughout the paper I consider only stellar/planetary magnetic fields that are in a 
"potential state." This terminology, which is common in the solar/stellar physics literature, 
refers to a magnetic field in which the current density J = everywhere within the stellar 
magnetosphere, and therefore the field B can be written in terms of the gradient of a 
magnetostatic scalar potential. 



II. STELLAR MAGNETIC FIELDS 

Surface magnetic field maps derived using the technique of Zeeman-Doppler imagin^^ 
have been obtained for a large number of stars of differing age and mass. Zeeman-Doppler 
imaging is a tomographic technique that allows stellar magnetic maps to be constructed by 
monitoring distortions in spectral lines recorded in circularly polarized light. Large scale stel- 
lar magnetospheres show varying degrees of complexity, from the simple well ordered almost 
dipolar field of a low mass star^^ (V374Peg, ~ 0.3 Mq) to the highly complex multipolar 
field of a high mass star^'^ (r Sco, ~ we show an example of a magnetic 

map derived using Zeeman-Doppler imaging of the forming star V2129 Oph (~1.35 Mq).'!^ 
Forming stars are still contracting under gravity. At this early phase of stellar evolution, 
the temperature/pressure in their cores is not yet sufficient for the fusion of hydrogen into 
helium to begin. The large scale magnetic field of V2129 Oph was found to be dominantly 
octupolar.fl^ Also shown in Fig. |2]is a numerical extrapolation showing the three-dimensional 
field topology constructed from the magnetic map (see Ref. [T7]for a review of the numerical 
technique). 

Motivated in part by the availability of observational data, which for the first time is 
allowing variations in stellar magnetic field topologies with stellar age, mass, and rotation 
period to be probed, new theoretical models that incorporate multipolar magnetic fields have 
been develop ed.'^'^ The basic properties of multipolar magnetospheres have been explored for 
many years by geophysicists. In Ref. IH] an expression for the field lines of an axisymmetric 
(axial) multipole of arbitrary degree i is derived (see also Refs. [25ll26l and Ref. EZI for an 
extension to the field lines of sectorial multipoles). In the context of this paper the order i 
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of a magnetic multipole can be thought of as the number of polarity changes in the surface 
field between the north and south pole of the star along a line of constant longitude (see 
Fig. |3]). By "equation of the field lines" I mean an expression of the form r = r(9) which 
describes the path (shape) of the field lines in a spherical coordinate system. The simple 
form r = r{6) is valid for field lines within a particular meridional plane (planar field lines). 
In this paper I restrict our attention to axial multipoles, which generate planar field lines. 
For example, an axial dipole {i = 1) stellar magnetic field has spherical components 

Br = Bl'^^'^(^)\ose (1) 



Be = ^Bl'^"'^(^^ysme (2) 

where i?* is the stellar radius and Bl'^°^'^ is the strength of the dipole at the rotation pole 
of the star.13 The field components are sketched in Fig. |4j 

Throughout this paper r and 6 are standard spherical polar coordinates, as illustrated in 
Fig. |4| with 6 measured from the rotation pole of the star; ^ = corresponds to the stellar 
rotation pole, and 6 = 7i/2 corresponds to the equatorial plane. The origin of the coordinate 
system is at the center of the star, and thus the stellar surface corresponds to r = R^,. Only 
azimuthally symmetric multipoles, known as axial multipoles, are considered in this paper, 
that is, B^ = 0. 

The differential equation describing the path of the field lines i^^'^ 

Br _ Be_ , . 

dr ~ rdd' ^ ' 

which, upon substituting for B^ and Bq from Eqs. ([T]) and ^ and integrating, yields the 

well known result (see, for example, Ref. | 



r = constant x sin^ 6. (4) 

If a particular closed field line loop of the dipole reaches a maximum radial extent of rmax 
in the stellar equatorial plane, where 6 = ti /2 (see Fig. |4]), then the integration constant in 
Eq. (|4]) is equal to rmax- Thus for a dipole the equation of the particular closed field line is 

r = rmaxsin^6'. (5) 

Different values of rmax correspond to different field lines. Closed field lines with larger 
values of rmax have footpoints at higher latitude (that is, at smaller co-latitude 6). As shown 



in Fig. [T] the field line footpoints are where the closed loop under consideration connects to 
the stellar surface. (Open field lines have only a single footpoint on the star.) 

Equation (|5| also describes the field lines of a point electric dipole, and Eq. ([s]) can be 
written in other coordinates for both the magnetostatid^^'^and electrostatic cases. The later 
expression is found in many electromagnetism textbooks such as Ref . EH their equation (2) in 
their section 1.08. The equation of the field lines for an arbitrary order finite electric linear 
multipole, for example, can be derived from the differential equation E^/dx = Ey/dy = 
EjdzW^ 

III. THE MAGNETIC FIELD COMPONENTS OF AXIAL MULTIPOLES 

The large scale magnetic field external to a star or a planet is generated due to dynamo 
action and the distribution of current sources within its interior. The region external to the 



star can be assumed to be source free, an approximation which is discussed in Sec. Ill A 
The spherical magnetic field components of stellar magnetospheres were derived in Ref. [17] 
by an expansion of the magnetostatic scalar potential using this assumption and the work 
of Refs. |2HI and |5] I do not reproduce the lengthly derivations here. 

In this paper I consider axial multipoles where I have chosen the space-fixed 2-axis of the 
star (the stellar rotation axis) as the symmetry axis of the multipole being considered (see 
Ref. [TTjfor details). As shown in Fig. |4]the field lines of axial multipoles are reflectionally 
symmetric in the x-axis, which lies in the stellar equatorial plane, and rotationally symmetric 
about the z-axis. The spherical field components of an axial multipole [B^ = 0) of order i 
are 

PiicosO) (6) 

where Bl'^°^'^ is the field strength of the multipole being considered at the stellar rotation 
pole, P* is the stellar radius, r is a point external to the star, P£i(cos^^) and Pi{cos9) are 
the m = 1 ith associated Legendre function and the ith Legendre polynomial respectively; 
i represents the multipole being considered, for example, £ = 1 is the dipole and i = 2 the 
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quadrupole. I define the associated Legendre functions as 



JTTt 

Pem{x) = {l-x'r/'—Pe{x), (8) 
wliere tlie Legendre polynomials are 

and X = cos 9. The associated Legendre functions are defined without the Condon- Shortley 
phase (a term of the form (—1)"^), as in Ref. [T71 The basic properties of the associated 
Legendre functions and the Legendre polynomials can be found in many texts, for example 
Ref. 123 Their functional form as pertains to the spherical magnetic field components is 
discussed in the appendix.!^ 

A. The effects of stellar outflows 

Equations ^ and ([T]) do not take into account the influence of outflows on the structure of 
large scale stellar magnetospheres. At some distance above the stellar surface, the coronal 
plasma distorts and pulls open the closed field line loops to produce a stellar wind. To 
account for this transition from the closed coronal field to the open interplanetary magnetic 
field for the Sun, Refs. [5] and [7] introduced the potential field source surface model. A 
potential field is one which is current free, that is, the current density J = everywhere 
within the stellar magnetosphere, and the field B can be written in terms of the gradient of 



a magnetostatic scalar potential that is, B = — V\l' (see Sec. IIIB). The source surface 
is a spherical surface of radius Rg at which the magnetic field is purely radial {Bg{Rs) = 
B^{Rs) = 0). This equipotential surface mimics the effect of the coronal plasma, which 
pulls open the large scale closed field lines (although the field may become radial inside this 
radius). In-situ satellite observations of the large scale solar heliospheric field within the 
region of space around the Sun where the solar wind has blown away the gas and dust that 
pervades the rest of interstellar space, suggest that the value Rs ~ 2.5 R© is reasonable for 
the solar corona. 

Several papers have examined the validity of potential field source surface model, and 
have generalized it by including current sheets outside of Rg and by considering a non- 
spherical source surface. The distortion of the field by the coronal plasma at the source 
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surface induces a current and invalidates the assumption that the large scale magnetosphere 
can be described as current free. Reference E2I found that the structure of the field close 
to the source surface obtained from magnetohydrodynamic (MHD) models that considered 
the influence of the coronal plasma on the field line shape could be reproduced by adop- 
tion of a prolate spheroid source surface (with major axis aligned with the solar rotation 
axis). The current-free assumption also often fails at reproducing the field topology of small 
scale field regions commonly detected on the Sun. However, for stellar observations, such 
small scale field regions remain well below the resolution achievable from the equivalent 
solar observations, and the current-free source surface model has been used extensively to 
produce three-dimensional models of the large-scale solar corona and stellar coronae via 
field extrapolation from observationally derived surface magnetic maps.^ ^^ * ^^ ' Reference |33] 
found that the large-scale field structure of the solar corona is adequately reproduced by 
the current-free source surface model in the sense that it produces magnetic field topolo- 
gies that are similar to those obtained from more complex and computationally intensive 
MHD field extrapolations. Reference [M] also found from their MHD simulations that the 
global structure of stellar magnetospheres show little departure from a potential state when 
evolved over several stellar rotations. The source surface boundary condition has thus been 
successfully applied to construct models of the magnetospheres of stars of various ages and 
spectral types,^ ° * ^^ * ^^ l as well as planetary magnetospheres.!^ Detailed comparisons between 
potential field source surface models and MHD models of global magnetic field topologies 
can be found in Refs. [33] and [T71 



B. The magnetic field components with a source surface 

For a spherical source surface of radius Rg the field components, Eqs. ^ and ([T]), are 



modified [see Eqs. ( 18 ) and ( 19 )]. In this section Laplace's equation is solved using separation 
of variables in spherical polar coordinates in the region R^, < r < Rg. From the solution I 
derive expressions for the magnetic field components of a multipole of order £ subject to the 
following conditions: currents within the region external to the star, r > R^,, and interior 
to the source surface, r < Rs, can be neglected, and therefore the field B can be expressed 
as B = — V\l/; the potential is azimuthally symmetric does not depend on 0, and hence 
= everywhere); the radial field component is specified at r = R^; and the polar field 
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component is zero at the source surface Rg, that is, Bg{Rs) = 0. 

Details of the derivation of the spherical field components with the source surface bound- 
ary condition are given in Ref. IT7| and only a brief summary is provided here. The large 
scale stellar magnetosphere is assumed to be static. Maxwell's equations for a magnetostatic 
field are 

V ■ B = (10) 
VxB = ^3, (11) 



where, assuming that currents can be neglected in the region external to the star, Eq. (11) 
reduces to V x B = 0. Hence we obtain Laplace's equation, V^\l/ = 0, whose general 
solution in spherical coordinates can be found in most texts on electromagnetism (see, for 
example, Ref. ESI), 

^ = 5^ 5^ [aen^r' + h^r-^'+''^] P,„(cos^)e^'"'^ {R. < r < R,). (12) 

£ m 

The coefficients and depend on the boundary conditions. Only the axial multipoles 
(with m = and hence oc d'^ /dcj) oc m = 0) are considered in this paper, for which a^o = 
a^, hiQ = bi and from Eq. ([s]) P^o(cos^^) = P^(cos6'). One boundary condition is to specify the 
radial field component at the stellar surface. For numerical field extrapolation models this 
condition comes directly from the observationally derived magnetic surface maps.'^'^'^ The 
radial component of the field across the stellar surface for a particular multipole of order i 
is given by Eq. ^ with r = R^, 

Br.{R,) = Bi^^°''Pe{cos9). (13) 

The other boundary condition is that the field at the source surface, r = Rg, becomes purely 
radial, that is, 

Be{Rs) = 0. (14) 
The field components for a particular axial multipole of order £ can be found from B = — V^' 



with Eq. (12) 



Be = - [a,/-' + 6,r-(^+2)] ^P^cos^). (16) 
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From Eq. (16) the boundary condition (14) is satisfied if 



be = -a,Rf+\ (17) 
which holds for any value of £.11331] jf substitute Eq. ( 17 ) into Eq. ( 15 ) , and apply Eq. ( 13 ) , 



we can derive expression for in terms of R^, and Rg. We then substitute this expression 



and Eq. (17) into Eqs. (15) and (16) and obtain general expressions for and Bq for the 



large scale magnetosphere with a source surface, 



B, 











f R* 




\ r 



1+2 



Pfi COS ( 



er 



2£+l 



{i + l)Rf+' 



1+2 



Pii (cos ( 



iRf+^ + + l)Rf+^ 
-(£ + l)r2^+i + (£ + l)i?f+i 



^R 



.21+1 



+ (£ + l)Rf+^ 



(18) 
(19) 



for R^ < r < Rg. To derive Eq. (19) we used the fact that P^i(cos^^) = — dP£(cos^)/d^ 
(see Eq. ([s]) with m = 1 and x = cos 9). The field components describe the structure of the 
magnetic field in the region exterior to the star/planet (r > i?*) and interior to the source 
surface (r < Rg). At radii beyond Rg, as for the Sun, the field is usually described by a 
Parker spiraP^ (the large scale solar field being dragged into a spiral shape by the outflowing 
wind), which is not considered in this paper. As Rg is decreased, more of the closed field 



line loops are converted to open field lines, and in the limit Rg — )■ oo, Eqs. (18) and (19) 
reduce to Eqs. ^ and ([T]). From the modified field components it is possible to derive the 
equation of the field lines. 



IV. EQUATION OF THE FIELD LINES 

The differential equation for the field lines of an arbitrary magnetic multipole i i^SEHSZl 

Bj. Bq Bfj) B . . 

dr rdO r sin Odd) ds ' 



where s is a coordinate along the path of the particular field line under consideration, and 
-B = |B| = (B"^ + Bg + B'^y^'^ is the magnitude of the field. For axial multipoles the equation 
describing the path of the field lines can be determined by solving the first two terms of 
Eq. (20).'^ Likewise, the equation for the field lines of a point electric linear multipole can 



be found in an analogous manner, see for example Appendix A of Ref. UHl either in spherical 
or Cartesian coordinates .'2322 
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We substitute Eqs. (18) and (19) into Eq. (20) and rearrange terms to obtain 

"-(£+l)r2^+i + (£ + l)p2m" 



{i + l)PKcos 
Pfl(cos6'' 



dr 
r 



£r2^+i + (£+l)p2m 



(21) 



To obtain the equation of a particular field line that has a maximum radial extent of 



(see Fig. |4j) Eq. (21) must be integrated from Tmax to some point r along the loop where 
R* < r < Tmax and P* < rmax < Rs- For closed field lines at r = rmax, Pr = 0, and the field 
only has a Bg component (see Fig. |4]). For open field lines rmax = 



Rg. The right-hand side 



of Eq. (21) is integrated by writing the integral as 



±1) 

+ r 



In {ir^'+^ + (3) 



+ 



f3r dr 



r2^+i(£r2^+i + 13)' 



(22) 



where /3 = (£ + 1)P^^"'"^. The inclusion of the additional r"^^ term in the second integral 
allows it to be simplified by the substitution a = ir'^^^^ + /3. The second term can be further 
simplified by partial fraction decomposition and directly integrated to give 



In 



^r^+x^ + + mf^' 

1^21+1 + + l)i?2£+l 



(23) 



The left-hand side of Eq. (21) is more complicated to integrate. We want to integrate 
from the maximum radial extent of a given field line to some point along the field line at r. 



I assume that at rmax, = Q, and the left-hand side of Eq. (21) is 

{i + l)Pi{cose) 
iQ Pii{cos6 



f 

Je 



-de. 



(24) 



Equation (24) can be integrated following the arguments in Ref. |2ll However, given the dif- 
fering definitions of the associated Legendre functions, and the fact that I also wish to eval- 
uate the integration constant, it is instructive to briefly summarize the method here. From 
Ref. ISHl Eq. (15), the Legendre polynomials are related to the P£i(cos6') and P(£_i)i(cos6') 
associated Legendre functions as 



P/(cos( 



cos 6 



Pn(cos( 



1 



'sin^ 



P(^_i)i(cos6'). 



Using this relation Eq. (24) becomes 

+ 1\ cos6' 
sin 6 



e 



1\ 1 P(^_i)i(cOSI 



sin 6* P£i(cos 



de. 



(25) 



(26) 



The first term can be integrated directly, but the second term cannot. From Ref.|38l Eq. (19), 

,dPei(cos 9 



— sinO- 



de 



-£ cos 6'Pa (cos 6') + (£ + l)P(£_i)i(cosi 
10 



(27) 



(28) 



I also note that 

[sm2^Pa(cose)] = 2 cos OPn (cos 6) + sin 6 "^^^^^^"^ ' 

By combining Eqs. (27) and (28) an expression for (£ + l)P(f_i)i(cos 6^) can be obtained, 
which, upon substitution into Eq. (26), and some manipulations. 



1 d 
sin 9 dO 



yields 



/e 



Icos^ l-|[sin2 0P,i(cos( 



de 



sin 9 i sin^ OPn (cos ( 



de. 



(29) 



Equation (29) can now be integrated directly, 

sin6'P£i(cosi 



In 



(30) 



sin 6P£i(cos 6) 

The modulus signs can be removed as the signs of P£i(cos6'), P£i(cosO), sin 6*, and sinG 



are always the same (see the appendix). Finally, by combining Eqs. (23) and (30), the 
solution of Eq. ( [2l| ) which is the equation of the field lines of an axial multipole with a 
source surface is found to be 



sin ^P^i(cos I 



lr2e+i + + l)p2£+i (£^2^+1 + + l)Pf +1) sin ePfl(cos 6) ' 



(31) 



where f = r/R^, Rg = Rg/R* and fmax = '^'maxZ-R*- Although Eq. (31) may lack simplicity, 
once a source surface location and the maximum radius of the field line of interest has been 

Different 



chosen, the right-hand side of Eq. (31) simplifies to constant x sin^P^i(cos 
field lines correspond to different values of this constant. Each region of closed field lines 
(see Fig. [T]) must be considered separately because each region of closed loops reach their 
maximum radial extent at a different value of G. For closed field line loops is the value 
of 6 where Pe{cosd) = (and consequently P^ = (see the appendix). For open field lines 
'"max = Rs, and G is the polar coordinate where a particular open field line crosses the source 
surface. 

The footpoints of a particular closed field line loop, that is, the co-latitudes 6^: where the 
closed loop attaches to the stellar surface (see Fig. [T]), are obtained by integrating Eq. (21) 

from Tmax 



to P*, which is equivalent to replacing r by 1 in Eq. (31) 



sin 6'*P£i(cos 6^,) 



V i+{i+l)Rf+^ 



sin 9Pfi(cos 9). 



(32) 
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There are two solutions 9^ of Eq. (32) corresponding to the two footpoints of the particular 
closed field line loop under consideration. 

The location of the source surface Rg controls the amount of open field relative to the 
closed field. Of particular interest is the amount of open flux through the stellar surface. 



By integrating Eq. (21) from rmax = -Rs to r = i?* [achieved by substitution of fmax = Rs 
and f = 1 into Eq. (31)] it is possible to determine the location of the boundaries between 
regions of open and closed field lines at the stellar surface which is where the largest closed 
field line loop connects to the star (for example, see the long-dashed field line in Fig. [T]), 

(2f 11^^+^ 

sm9^Pn{cos9^) = sin eP,i(cos 6). (33) 



From Eq. (33) it is possible to derive an analytic expression for the open flux through the 
stellar surface, which is an important quantity for models of stellar rotational evolution. Due 
to angular momentum loss via outflows, the rotation rates of low mass stars slow during the 
hydrogen burning phase of their evolution.'^Sl The angular momentum that can be removed 
by a stellar wind depends on the amount of open flux. 

A. Field line plots 



The equation of the field lines, Eq. (31), makes possible a visualization of the magnetic 
field topologies of arbitrary magnetic multipoles of order i, which incorporate regions of 
both closed and open field lines. Figure [5] shows the field lines of the three lowest order 
axial multipoles, the dipole {£ = 1), the quadrupole (£ = 2) and the octupole (£ = 3). 
For a multipole of order i there are 2i regions of closed field line loops around the entire 
circumference of the star in each meridional plane (planes with = constant). Additional 
properties of axial multipoles are discussed in the appendix. The path of the field lines 



illustrated in Fig. 5 are obtained by solving Eq. (31) as follows. Suppose we are interested 



in a particular closed field line loop that has footpoints 9i and 92, and which reaches its 



maximum radial extent at G. The footpoints 9i and 92 are obtained by solving Eq. (32) 



For a given region of closed field lines for an axial multipole of order i and given values of Rs 



and fjnax, the right-hand side of Eq. (32 ) is a constant, say rj. The footpoints of the closed field 
line 9^: can then be found by finding the roots of the function /(0*) = sin ^^*P; i(cos 0*) — 
rj using a standard root finding algorithm. There are two roots 9^, corresponding to the 



12 



footpoints 6i and 62- Because axial multipoles are rotationally symmetric about the 2;-axis 
and reflectionally symmetric in the x-y plane, I only consider values of 6 between zero and 
7t/2. Odd i multipoles have a region of closed field lines that span the stellar equatorial 



plane, G = 7r/2 (see Fig. pi), and for this region there is only one solution 6^, of Eq. (32) for 



< 6 < 7t/2. Once the di and 62 of a closed loop have been determined, it is straightforward 



to solve for the path of the field line by solving Eq. (31 ), again using a root finding technique 
to find f at each value of 6 between 61 and 92- Equation (31) can also be used to trace the 
paths of open field lines. For open field lines, fmax is replaced by Rg in Eq. (31 ), and 9 is the 
polar angle where a particular open field line arrives at the source surface. For example, for 
a particular multipole that has closed field lines which reach their maximum radial extent 
at 6i, where i = 1,2, . . ., the open field lines arrive at the source surface (r = Rs) at all 
angles between and 7i/2 except Gj. 

V. SUMMARY 

Over the past five years the current generation of optical spectropolarimeters has allowed 
stellar magnetic field topologies to be probed in unprecedented detail. Large observational 
programs have provided maps of the magnetic fields of stars spanning a range of masses, 
rotation rates, and evolutionary stages.'^ In particular, magnetic maps of forming solar- 
like stars can now be obtained.^ The large scale field topologies vary from simple dipoles 
to more complex magnetic fields consisting of several multipole components.'^ Given the 
availability of the observational datasets, new models that incorporate magnetic fields are 
under development such as models that consider how the magnetospheres of forming solar- 
like stars interact with their surrounding planet-forming fligl^s.liillllMI 

I have discussed how the solution of Laplace's equation using separation of variables can 
be applied to construct models of stellar magnetospheres where the influence of outflows on 
the global field structure is considered. With the assumption that currents have a negligible 
effect on the large scale structure of stellar magnetic fields in the region between the stellar 
surface and the source surface, the field components have been derived by solving the Laplace 
boundary value problem, thus by passing the more complicated problem of solving Poisson's 
equation. Basic properties of multipole magnetic fields have been discussed, and from the 



field components, Eqs. (18) and (19), an analytic equation for the field lines of an arbitrary 
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axial multipole with a source surface has been derived. The resulting expression, Eq. (31), is 
straightforward to solve using a root finding algorithm and its use only requires knowledge 
of the roots of the m = 1 ith associated Legendre function P^i(cos^). 

I have concentrated on deriving the equation of the field lines for a modified magnetic 



multipole in spherical coordinates by solving Eq. (20). Spherical coordinates are the most 
natural coordinate system to use in the description of stars and planets and their imme- 
diate environments. Other authors have derived analytic expressions for the field lines of 
axisymmetric multipoles without the source surface boundary condition in spherical polar 
coordinates.'^^mS] rpj^g external field line equation for the magnetic multipoles is identical to 
that obtained for the field lines of point electric linear multipoles. (However, the internal 
and contact field lines of magnetic and electric multipoles of the same order differ.'^ For 
example, the path of the electrostatic field lines for a linear quadrupole in spherical coor- 
dinates, can be found in Appendix A of Ref. [181 For the electrostatic case it is much more 
common for the path of the field lines to be determined in Cartesian coordinates. Ref- 
erence [21] provides a Cartesian expression for the field lines of an arbitrary linear electric 
multipole. Several other overviews describe algorithms for simple computer programs that 
can be adapted to visualize magnetostatic or electrostatic fields.'^^'^ZEH] 

The magnetospheric structures considered in this paper have certain limitations. The 
fields are assumed to be axisymmetric. The assumption of symmetry allows the analytic 
equation for the field lines to be derived. In reality, stellar magnetic fields can be highly 
complex, with many high order field components. (See Ref. |T2I for a review of recent re- 
sults on stellar magnetic topologies.) A different approach is to model stellar magnetic fields 
numerically via extrapolation from observationally derived magnetic maps. '^'^'^ Both the nu- 
merical approach and the analytic work in this paper assume that the stellar magnetic fields 
are current free. As discussed in Sec. Ill A , the global field topologies obtained in potential 
field source surface models approximately match those obtained from more computationally 
complex magnetohydrodynamic mo dels .'^^'^ 

In this paper I have assumed that the ith order multipole moment symmetry axis is 
aligned with the stellar rotation axis, and that both lie in the same stellar meridional plane 



(planes with = constant; for example, the x-z plane). However, Eqs (31)-(33) apply gener- 
ally to tilted multipole symmetry axes (that is, large scale magnetospheres tilted arbitrarily 
in the polar and azimuthal directions with respect to the stellar rotation axis, assumed to 
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be the z axis) with appropriate coordinate and vector frame transformations. 

The equation for calculating the co-latitude of the footpoints of the largest closed field line 
loop within each region of closed field lines can be used as a basis for calculating the amount 
of unsigned open (or closed) flux relative to the total flux (open plus closed) through the 
stellar surface. Such detailed calculations, from which models of stellar rotational evolution 
can be developed (for example, Refs. |lT]and|l2]), are deferred to a future paper. 



APPENDIX: BASIC PROPERTIES OF AXIAL MULTIPOLES 



Consider a multipole of order i with a single closed fleld fine extending to a radius rmax 
above the stellar surface within each region of closed fleld lines. As Q increases from the 
rotation pole of the star {Q = 0) to the equatorial plane {6 = 7r/2), the sign of and Bg will 
vary according to whether P£{cos6) and P£i(cos^) (respectively) are positive or negative. 
At the rotation pole of the star Pi{cos6) = 1 and P^i(cos6') oc sin 6' = 0. The fleld therefore 
only has a radial component, and Bg = 0. As we move from the rotation pole toward the 
highest latitude closed fleld line footpoint, both P^(cos6') and P^i(cos6'), and therefore Br 
and Bg, are positive. As 6 is further increased, P£(cos6'), and therefore Br, changes sign 
every time we pass through the values of 6 (B in our notation) where the closed fleld line 
loops reach their maximum radial extent (for odd i multipoles there is a region of closed 
fleld lines that reach there maximum radial extent in the stellar equatorial plane; for this 
region 9 = tt/2). Pgi{cos9), and therefore Bg, changes sign at values of 9 within the open 
fleld line regions that occur between the regions of closed fleld lines. Hence, for a particular 
closed fleld line with footpoints at 6i and 62 and which reaches its maximum radial extent 
'"max at O, P^i(cos^^) has the same sign at each point along the loop. For this reason the 



modulus signs in Eq. (30) can be removed. 

These ideas are best illustrated by considering an example of a particular multipole. 
Figure [6] shows a hexadecapole (i = 4) magnetic fleld. For clarity, only one loop within 
each region of closed fleld lines is drawn. The high latitude closed loop is assumed to have 
footpoints at 61 and 62, and the lower latitude closed loop has footpoints at 63 and 64^. 
The closed fleld lines in the upper region reach their maximum radial extent at Gi, and 
the closed fleld lines in the lower region do so at O2. The closed loops are not symmetric 
about the Br = lines because the Legendre polynomials are not symmetric about their 
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roots. This property means that the angular difference between each footpoint and where 
the particular region of closed field lines reach their maximum radial extent is not the same; 
for the particular hexadecapole field line Qi — 9i ^ 62 — Qi- The exception to this is the 
equatorial region of closed loops for odd ^ multipoles, which are symmetric about the stellar 
equatorial plane. 
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FIG. 1: (color online) A star of radius i?*, shaded gray, with an axial dipole magnetic field with 
a source surface (the large dashed circle) of radius Rs- The source surface allows for the inclusion 
of regions of open field lines, along which a stellar wind is launched, in addition to the regions 
of closed field lines. Each closed field region is surrounded by regions of open field lines. The 
field lines are shown in light gray. The solid black lines illustrate particular field lines, with their 
footpoints on the stellar surface denoted by the solid black circles. Open field lines have a single 
footpoint, and closed field lines have two footpoints on the star. The field lines within each region 
of closed field lines reach their maximum radial extent of rmax along a line of constant polar angle. 
Different values of rmax correspond to different field lines within each closed field region. The 
long-dashed line denotes the closed field line loop with rmax = Rs- The footpoints of this closed 
field line represent the boundary at the stellar surface between regions of closed and open field. 
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FIG. 2: (color online) (a) A flattened polar projection showing the radial field component of the 
surface magnetic field of the forming star V2129 Oph (~ 1.35 M©).'^^! The stellar equator is shown 
as the bold circle, with lines of constant latitude separated by 30° as the dashed lines. There is a 
positive field spot slightly offset from the pole, surrounded by a ring of negative field, which itself is 
surrounded by a ring of positive field below the equator. Fluxes given in Gauss. Numbers and tick 
marks around the circumference denote the rotation phase and phases of observation respectively. 
The large scale magnetic field is found to be dominantly octupolar, which is apparent from the 
(b) extrapolation of the three-dimensional coronal magnetic field, which is constructed from the 
magnetic map. Closed field lines are shown in white, with open field lines, along which a stellar 
wind could be launched, shown in blue (see the online version for colors). The large scale field 
approximately resembles a tilted dipole, and the medium scale field resembles an octupole with 
three rings of closed field. Details of the observational techniques used to construct stellar magnetic 
maps and the numerical field extrapolation model can be found in Refs. [T^] and [T7\ respectively. 
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FIG. 3: (color online) The first three lowest order multipoles, the dipole {i = 1), the quadrupole 
(i = 2) and the octupole {£ = 3). The multipole moment symmetry axis, denoted by \ii for 
the dipole, for the quadrupole, and yi^ for octupole, is assumed to be aligned with the stellar 
rotation axis, denoted by ft. The plus/minus signs denote regions of positive/negative field. The 
order of a multipole i is the number of polarity changes in the surface field between the north and 
south pole of the star along meridians (lines of constant longitude); i is the number of roots of 
Pi{cos9) between ^ = and 6 = tt. A multipole of order i has 2i regions of closed loops around 
the entire star in each meridional plane (planes with (p = constant). In three dimensions the closed 
field regions form rings of closed field around the star, with i rings for a multipole of order I. 
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FIG. 4: A field vector B decomposed into the radial Br and polar Bg components at a point 
along a field line a distance r from the center of the star (shaded gray) at a co-latitudc of 9. 
The field components are used to illustrate their definitions and are not to scale. The field line 
reaches a maximum radial extent of rmax- Only the first quadrant is shown because the magnetic 
fields considered in this paper are reflectionally symmetric in the horizontal axis and rotationally 
symmetric about the vertical axis. 
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FIG. 5: (a) The field lines of a dipole = 1), (b) a quadrupole {I = 2), and (c) an octupole 
[i = 3) with a source surface at Rg ~ 3.4 R* plotted as the dashed line. The star is shaded in gray. 



The shapes of the closed field lines are calculated for each multipole from Eq. (31). The shapes 
of the open field lines are calculated in a similar way from i?* to by considering the range 
of co-latitudes Q not occupied by the footpoints of the closed loops. For the quadrupole and the 
octupole (and other higher order multipoles) there are regions of open field at lower latitudes along 
which a stellar wind could be launched. The magnetic fields are rotationally symmetric about the 
z-axis and reflectionally symmetric in the x-y plane. 
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FIG. 6: (color online) A star (shaded gray) with an axial hexadecapole magnetic field {I = 4, 
m = 0). Q denotes where the closed field lines reach their maximum radial extent, and 9 indicates 
the locations of the footpoints of the closed loops. For clarity only two closed field line loops 
within each closed field region are shown (solid lines), and no open field lines are shown. The 
field is rotationally symmetric about the z-axis and reflectionally symmetric about the x-axis, 
and therefore we need consider 9 values only between zero and tt/2. The field is purely radial 
at the stellar rotation pole, in the equatorial plane, and along a line between the closed field 
lines shown here (and by symmetry in the other three quadrants). Because Br oc P£{cos9) and 
Bg oc Pgi (cos 9) , the values of 9 where Br and Bg are zero correspond to the roots of the particular 
Legendre polynomial and associated Legendre function, respectively. 
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